Nonlinear conductance in a ballistic Aharonov-Bohm ring 
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The nonlinear electronic transport properties of a ballistic Aharonov-Bohm ring are investigated. 
It is demonstrated how the electronic interaction breaks the phase rigidity in a two-probe mesoscopic 
device as the voltage bias is increased. The possibility of studying interference effects in the nonlinear 
regime is addressed. The occurrence of magnetic field symmetries in higher order conductance 
coefficients is analyzed. The results are compared with recent experimental data. 
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Introduction: The Onsager-Casimir [l[ reciprocity 
relations are ubiquitous in systems driven out-of- 
cquilibrium in the linear response regime. In classical 
systems, the quest for related symmetries in the non- 
linear regime is a subject of increasing interest with po- 
tential applications in numerous problems 0]. In meso- 
scopic physics, the Onsager symmetries enforce, for in- 
stance, the linear conductance in two-terminal devices 
to be an even function with respect to magnetic field 
inversion Often described as phase rigidity, this pre- 
vents the use of such devices for quantum interferometry 
0| . Recently, experimental groups have studied how the 
Onsager-Casimir relations are broken in nonlinear trans- 
port as the source-drain bias is increased in a variety 
of open systems, such as quantum dots 0, S 0|> carbon 
nanotubes [|| , and quantum rings 0, [l(| • A large bias in- 
duces a rearrangement of charges, which are subjected to 
the Coulomb interaction. Far from equilibrium, interac- 
tions break the Onsager-Casimir symmetries [HI [l2l. [l3| . 
Theoretical effort has been devoted to infer the effective 
electron-electron interaction parameters from the statis- 
tical analysis of nonlinear transport experiments [T3 |. 

A recent study of nonlinear transport in a two-terminal 
ballistic Aharonov-Bohm (AB) ring as a function of an 
applied perpendicular magnetic field B and a bias V [9| 
unveiled non-statistical puzzling symmetries. By expand- 
ing the current 1(B) in powers of V, 



G (1) V + G (2) V 2 + G (3) V 3 + 



(1) 



nonlinear conductance coefficients of different orders were 
analyzed. It was observed that the even conductance 
coefficients G' 2 ™' are neither even nor odd with respect 
to the magnetic field. In contrast, surprisingly, the odd 
ones, 

G (2n+1) 

are even functions of B raising the ques- 
tion whether this observation is an indication of a new 
fundamental non-equilibrium symmetry. Ref. also 
shows that phase rigidity is lifted in the nonlinear regime. 
The accumulation of an extra phase in one of the quan- 
tum ring's arms is mimicked by using lateral gates that 
locally modify the electronic density. The phase shift of 
the AB oscillations in G^ as a function of lateral gates 



voltages shows a smooth variation with a slope propor- 
tional to the arm length, which is a promising result to- 
wards using AB rings as two-terminal interferometers. 

In this paper we introduce a simple model that allows 
for a semi-analytic description of the nonlinear transport 
in a ballistic AB ring. We analyze the nonlinear con- 
ductance oscillation phase shifts and the symmetries of 
the nonlinear conductance coefficients, both issues ex- 
perimentally addressed in Ref. 0. We discuss the phase 
rigidity in the G^ conductance coefficient and present 
an explanation for the observed even parity in the G^ 
and G^ conductance coefficients. 

The method: Let us consider the electronic transport 
in a two-probe mesoscopic ring connected to reservoirs, 
a = 1,2, both at a temperature T. A voltage bias applied 
to the reservoirs drives the system out of equilibrium and 
causes a current flow. In the absence of inelastic pro- 



cesses, the electron current on the lead a reads [111. \14\ 



la = 



2 

E 

0=1 



dE fp(E)A a p(E;{U(r)}) , (2) 



where fp(E) = f(E - ft - eV ) = [ e (E-»a-eV )/k B T + 
1] , fcB is the Boltzmann constant, and Vg is the voltage 
applied to the reservoir connected to the lead (3 (mea- 
sured with respect to the equilibrium chemical potential 
fxo). Due to current conservation, I = I\ = — J2. 
The transmission coefficient A a p is given by [ll[ 



A af3 (E;{U(r)}) = Tr[l a 6 a (3 



(3) 



where S a p denotes the scattering matrix with lines (rows) 
associated with the channels a (b) at the contact a 
(/3), l a is the identity matrix, and the trace runs over 
all open channels in a and (3. The scattering matrix 
S a p(E, {U(r)}) is a function of the electron energy and 
a functional of the electrostatic potential U(r) that is 
established in the conductor due to the voltages {Vj}. 

It is convenient to expand AU(r) = U(r) — U eq (r) as 
a power series of the applied voltages, namely, 



a a,f3 



(4) 



2 



with the characteristic potentials u a p...(r) defined as 
( d 



u a p-{v) 



d 



\dV a dVp 



U(r) 



{y T }=o 



In line with Ref. @, we take Vi = V/2 and V 2 = -V/2. 

The linear conductance, G' 1 ) = dI/dV\v=o, is 
given by the Landauer formula, namely, — 
(2e 2 /h)JdE(-d E f )(A n - A 12 )(E;{U cq (r)})/2, where 
fo(E) — f(E — jLto) and U eq is the equilibrium electro- 
static potential for V = 0. 

Simple symmetry arguments show that G^(B) = 
G^ 1 -* (—£>). In contrast, there is no general principle 
that predicts reciprocity relations for the nonlinear con- 
ductance coefficients. The inspection of higher powers 
of V in Eq. (JTJ) gives any G^ in terms of the char- 
acteristic potentials and functional derivatives of the 
transmission coefficients with respect to U(r). While 
quantities like 8A a p/5U(r)\y = Q are expressed in terms 
of S a p(E; {C/ oq (r)}), the u a 0—'s encode the information 



about electronic many-body interactions 1J, [15| . 

The conductance coefficient G^ is written as G^ 
e 3 /2h dE (-d E f Q ) T« with [ll| 



= / dr [«i(r)-« a (r)] 



A, 



eSU(r) 



(6) 



At the simplest approximation level, u a (r) is obtained 
from a Hartree equation, whose source terms are the in- 
jected electron density and the corresponding induced 
charge in the conductor, the latter related to the Lind- 
hard function [l4j ] . Using the Thomas- Fermi approxima- 
tion for the polarization, one writes 11 , ll| 



- Au a (r) + 4ne 2 u a (r) 



dn°™ 



(r) 



dE 



dE 



(7) 



where the injectivity is given by 
d<(r) 



dE 



dE 



dfo 
dE 



/3 



Tr 



at ^ S /3a 

0a e5U(r) 



H.c. 



(8) 

and the emissivity, dn c ™(r)/dE, is obtained by exchang- 
ing a and f3 in S at Eq. ([5]) [14j. Injectivity and emis- 
sivities are related by symmetry: S a p(B) = [S a p(— B)] T 
readily gives that dn™(B)/dE = dn c ™(-B)/dE, often 
called microreversibility relation [Til ]. Gauge invariance 
requires that J2 a u a{ r ) — 1) leading to ^ Q dn^(r)/dE = 
E a ^ n ™( r )/^ = dn(r)/dE, the density of states. 

Since A a p and 5A a p/5U(Y)\v=o show the same sym- 
metries with respect to B, one concludes that the char- 
acteristic potentials are formally responsible for the vio- 
lation of the Onsager relations G' 2 ' . The same reasoning 
holds for higher conductance coefficients. 

Expressions for the coefficients G^ are known to ar- 



bitrary order [15|. G^ 



/ShJZ.dEi-dEfo)^), 



necessary for what follows, reads 
5(A n - A X2 ) 



(5) = dvu^(r) 



e5U(r) 



dr / dr' 



u (i)( r ) u (i)( r ') + i 



2 



6 2 (A 11 ~A 12 ) 



eSU(r)eSU(r') 



(9) 



v=o 



where = u\ — u 2 and = u\\ + u 22 — u\ 2 — u 2 \. 
uW is the solution of 0, HI] 

d 2 n 



A + 47re 2 ^ ) = -Aire 3 (u^) 2 + 1 

dill 



dE 2 



m ( d 2 nf d 2 nf\ , . 
- 2U( ) d^)^ 

which does not show any simple universal symmetry. We 
conclude that the symmetry experimentally observed in 
Q{n>i) are S y S tem specific, a subject we investigate next. 

Model: We address the electronic transport through 
an AB quantum ring using the single-channel model put 
forward in Ref. (l7| . The scattering problem is defined 
by the electron flow at the vicinity of the contacts as 
sketched in Fig. [T] Incoming electron wave function am- 
plitudes are indicated by unprimed Latin indices, while 
outgoing ones are denoted by primed letters. For in- 
stance, a a (a' a ) is the wave function amplitude of an in- 
coming (outgoing) electron at the contact a. 




FIG. 1: Two-terminal Aharonov-Bohm ring [17J: Arrows 
show the flow of the electron wave function amplitudes at 
the vicinity of the contacts, which are indicated by triangles. 

The scattering by the contacts connecting the leads to 
the ring is described by the scattering matrix S^ a \ As- 
suming that (i) the effect of the magnetic field on the 
electrons is negligible on the vertex length scale, (ii) the 
scattering amplitudes are symmetric with respect to the 
two branches of the ring, and (in) is real (the com- 
plex part of the S matrix can be introduced in the prop- 
agation at the arms of the ring); one writes [ijj 




/ i i * \ r/2 1/2 

1/2 



V 



,V2 



A a 



Xa 





( a a \ 


) 









(11) 



with Xa = -(Vl-2e Q - l)/2 and £ Q = -(^1 - 2e Q + 
l)/2. The parameter e a S [0,1/2] tunes the reflection 



3 



at the vertex a. The contact is closed for e a = and 
maximally open for e a — 1/2. When — the 
contacts are identical and the ring has, by construction, 
a reflection symmetry. 

The electron propagation is described by a transfer ma- 
trix which, for the upper arm, renders hi = e l< - kil ^^ 1 ' l b' 1 
and b' 2 = e^ kei+< ^ 1 'bi, where k is the electron wave vec- 
tor, l\ is the arm length, and <j>\ is the AB phase accumu- 
lated by the electron while traversing the ring from vertex 
1 to 2. A similar transfer matrix can be written for the 
amplitudes c and d by replacing l\ — > £ 2 and <f>i — > —02. 
The full scattering matrix S a p is readily obtained from 
and the above described transfer matrices. 

The electronic interactions are accounted for in 
SS a f}/8U(r) as follows. We place a 8- function scatterer 
at the position r along the ring, calculate the modified 
transfer matrices, and take the limit 

SS aP S aP [U(r') + n8(r' - r)] - S a0 [U(r')} 

SU(r) f™o n 

(12) 

8S a p/8U{r) is used to compute injectivities, emissivities, 
and 8A a p/8U{v). 

To relate our numerical findings to experiments, we 
choose the quantum ring diameter and set the effective 
electron mass to keep about 50 electrons below the Fermi 
energy. We present results for zero temperature. 

Results: Figures [2ja) and^b) illustrate typical injec- 
tivities calculated with our model. They present discon- 
tinuities at the contact positions, oscillations of period- 
icity about 1/kp around plateaus of different amplitudes 
for each arm. Oscillations in dn™(B)/dE are Friedel- 
like fringes due to trapped charges coming from contact 
a while the plateaus are related to the charge flow also 
coming from contact a. By computing the emissivity we 
verify microreversibility dn™ / 'dE(B) — dn c ™ /dE{—B). 
We also check that £ a dn™/dE{B) = J2 a dn c ™/dE{B). 

Since our primary goal is a qualitative understanding, 
we do not attempt to microscopically model the screening 
effects of a specific device. Instead, we use the contact 
interaction approximation and algebraically solve Eq. (JT]) 
to write the characteristic potentials as 



4ire 2 dn™/dE 
■? + 47re 2 dn/dE' 



(13) 



where n is a constant that characterizes the electron in- 
teraction strength. In the one-dimensional case consid- 
ered here, k is dimensionless. We set n = 0.1/ (kp£), 
with £ — i\ +£ 2 - Illustrative characteristic potentials are 
shown in Figs.[^c) and (d). 

Inserting u a (r) and 8A a p/8U(r) in Eq. ([6]) we obtain 
Qi 2 ), Figure [3fa) shows G^ versus magnetic flux (f>. 
Due to the system spatial symmetry introduced by taking 
= S^ 2 \ is odd in magnetic field. This peculiar 
nonlinear phase rigidity disappears as the spatial symme- 
try is broken. To realize this, we transform the scattering 
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FIG. 2: Injectivity at a = 1 (a) and a — 2 (b) for <j>/<j)o = 
±7r/4. The variable x G [0, 1] describes the position along 
the ring starting at the contact 1. Contact 2 is placed at 
x = 1/3. Characteristic potential, u^(x) for a symmetric 
bias and magnetic flux (c) <j>/<j>a = tt/4 and (d) 4>/(/>o = — 7r/4. 
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FIG. 3: Nonlinear conductance G' 2 ' as a function of <j>/4>o 
for both (a) a symmetric ring and (b) an asymmetric ring, 
(c) Nonlinear conductance for an asymmetric ring as a 
function of magnetic flux for different values of 8 = AMi. An 
offset was introduced to separate the curves, (d) Phase of the 
first harmonic of G^ as a function of the phase 8. 



matrix of a = 2 as -> R- 1 (6)S^ R^O), where R x (6) 
is the matrix representation of a classical rotation by 9 
around the x-axis. We observe that as increases the 
symmetry of G^ is gradually broken. Fig. [3{b) shows 
G< 2) for 9 = it/6. 

Phase locking in two-terminal devices is lifted as the 
source-drain bias is increased. The phase shift in the AB 
oscillations of G^ were experimentally studied by con- 
trolling the electron density at a selected arm [9|. Wc 
model this by changing the wave vector at one of quan- 
tum ring arms, namely, k — > k + Ak. This is equivalent to 
increasing the accumulated phase of the electrons flowing 
through this arm by 8 — Ak£i. Figure [3Jc) shows the 
conductance G^ for several values of 8. 



4 



In experiments, dephasing suppresses the higher har- 
monic contributions to that correspond to electron 
trajectories that go around the ring more than once. 
Since our treatment does not include decoherence effects 
we extract the first harmonic, G , of G^ to interpret 
the experimental data. Using the formula tan(A^) = 
(G( 2 )(0)sin(0/0 o ))/<G (2) (^)cos(^/(/.o)) @ we evaluate 
Acf> for different values of 6. Figure 03d) collects the val- 
ues corresponding to the curves shown in Fig. [3jc) . The 

(2) 

dependence of G on S and the magnetic flux (f>/(f>o is 
shown in Fig. 2J 
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FIG. 4: [Color online] G (2) , the first harmonic of G (2) , as a 
function of magnetic flux (j)/4>o and S. 

In Fig. El^d) we observe a smooth monotonic depen- 
dence between Acfi and <5. Although the dependence is 
not linear, except for small 6, it is in accordance with 
the experiment Q. However, since the perturbation 5 
modifies both, the characteristic potential and the func- 
tional derivative of A a p in a non-trivial way, a monotonic 
behavior is not generic. Varying the model parameters, 
we also find cases where A<p shows fast jumps as well as 
non-monotonic behavior as a function of 5. We conclude 
that despite eliminating the phase rigidity by applying a 
large bias, it is still a difficult task to use the nonlinear 
transport regime for interferometry studies. 

We now turn our attention to the parity of G^ n ' with 
respect to B. We already concluded that, in general, G^ 
is neither even nor odd in B. Let us present a system spe- 
cific scenario that explain the experimental observations 
Q: Most inelastic processes are very weakly dependent 
on the magnetic field. For a large bias, those give even 
parity contributions to the I — V characteristics. Hence, 
we focus on the odd parity contributions. When the con- 
tacts are equivalent, the spatial symmetry imposes that 
ui(B,r) = U2{— B,r) and the quantum coherent part of 
G^ 2 ) is odd in B. Similar arguments, applied to in 
Eq. (fTU)) . lead to a purely even G' 3 L In this scenario, 
inelastic processes even in B explain why G^ 2 ' does not 
have a definite parity and G^ has. The fact that the ring 
is quite open, favors the hypothesis of nearly equivalent 
contacts. Further support for this picture comes from the 
analysis of the h/2e oscillations observed at low magnetic 
fields [9(. Those can be attributed to paths going twice 



around the ring or to time-reversed paths each travers- 
ing the ring only once [lij]- While the first ones have no 
define B parity, the latter cause the Altshuler-Aronov- 
Spivak oscillations, which are usually more robust and 
by nature even functions in B. Hence, the absence of 
AAS oscillations in even part of the measured G^ and 
G^ is consistent with the hypotheses of a nearly sym- 
metric ring. 

In summary, we presented a theoretical approach that 
allows for a semi-quantitative discussion of nonlinear 
transport properties of ballistic mesoscopic devices, in 
particular of AB rings. We described how phase rigidity 
is broken in two-terminal devices as the voltage bias is in- 
creased. We also found that the experimentally observed 
symmetry G^ 2n+1 \B) = G (2n+ls >(-B) is not generic and 
can be explained in terms of systems specific features. 
Further experiments, exploring other control handles and 
geometries, can greatly contribute to settle this problem. 

We thank M. Biittiker, K. Ensslin, and R. Leturcq for 
valuable discussions and acknowledge financial support 
from the Brazilian funding agencies FAPERJ and CNPq. 

Note added: Related phenomena were very recently 
reported on a different system: An AB-ring with dis- 
crete level spectrum due to an embedded quantum dot 
(QD)[l^]. Here, the weak tunneling strengths make the 
asymmetry in the u's much weaker than in our paper and 
many-body effects in the QD can play a dominant role. 
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